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1 ATTORNEY DOCKET NO. 84432. 
2 

3 A METHOD FOR ESTIMATING THE PROPERTIES OF A 

4 SOLID MATERIAL SUBJECTED TO COMPRESSIONAL FORCES 

5 • 

6 STATEMENT OF GOVERNMENT INTEREST 

7 The invention described herein may be manufactured and used 

8 by or for the Government of the United States of America for 

9 governmental purposes without the payment of any royalties 
10 thereon or therefore. 

11 

12 BACKGROUND OF THE INVENTION 

13 (1) Field of the Invention 

14 The present invention relates to a method to measure (or 



15 estimate) the complex frequency- dependent dilatational and shear 

16 wavenumbers of a single slab of material subjected to large 

17 static compressional forces. More particularly, this invention 

18 provides a method to determine complex dilatational wavespeed, 

19 complex shear wavespeed, complex Lame constants, complex Young's 



20 modulus, complex shear modulus, and complex Poisson's ratio. 

21 (2) Description of the Prior Art 

22 Measuring the mechanical properties of slab-shaped materials 

23 are important because these parameters significantly contribute 

24 to the static and dynamic response of structures built with such 

25 materials. One characteristic that most elastomeric solids 
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1 possess is that, when they are subjected to large static forces 

2 (or pressure), their rigidity changes. Materials that have one 

3 set of mechanical properties at a pressure of one atmosphere can 

4 have very different properties when subjected to increased 

5 pressure. The ability to determine the pressure dependence of 

6 material properties is extremely important for modeling the 

7 behavior of systems comprised of these materials. 

8 . Resonant techniques have been used to identify and measure 

9 longitudinal and shear properties for many years. These methods 

10 are based on comparing measured eigenvalues to modeled 

11 eigenvalues and calculating the resulting material properties. 

12 These methods do not account for static pressure or large 

13 compressive forces. Additionally, they typically require long, 

14 slender materials to perform the measurement process. Comparison 

15 of analytical models to measured frequency response functions are 

16 also used to estimate stiffness and loss parameters of a 

17 structure. When the analytical model agrees with one or more 

18 frequency response functions, the parameters used to calculate 

19 the analytical model are considered accurate. If the analytical 

20 model is formulated using a numerical method, a comparison of the 

21 model to the data can be difficult due to dispersion properties 

22 of the materials. These methods do not take into account large 

23 compressive forces. 

24 In the prior art, some efforts have been made to measure 

25 material properties. under large pressures. These methods consist ■ 
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of placing materials in pressurized settings, insonifying them, 
and then measuring their response. These methods are difficult 
because they have to be conducted under great atmospheric 
pressure that can adversely affect the instrumentation. Safety 
issues can also arise in connection with laboratory testing at 
extreme pressures. Finally, a mass loaded long thin rod has been 
studied with respect to the bar wavespeed and corresponding 
Young's modulus; however, this work does not investigate shear 
motion. 

Accordingly, there is a need for a method of measuring 
mechanical properties of slab-shaped materials placed under 
pressure. 

SUMMARY OF THE INVENTION 
Accordingly, in this invention, a method to measure the 
complex frequency- dependent dilatational and shear wavenumbers of 
a material under a static compressional force is provided. The 
material is first vibrated in both vertical and horizontal 
directions while obtaining transfer functions in each direction. 
The two transfer functions are combined with a theoretical model 
to estimate a dilatational wavenumber and a shear wavenumber. 
The wavenumbers can be combined to give the complex dilatational 
wavespeed, complex shear wavespeed, complex Lame constants, 
complex Young's modulus, complex shear modulus, and complex 
Poisson's ratio. 



1 BRIEF DESCRIPTION OF THE DRAWINGS 

2 A more complete understanding of the invention and many of 

3 the attendant advantages thereto will be readily appreciated as 

4 the same becomes better understood by reference to the following 

5 detailed description when considered in conjunction with the 

6 accompanying drawings wherein: 



7 FIG. 1 shows apparatus for measurement of transfer functions 

8 in a vertical direction according to the current invention; 

9 FIG. 2 shows apparatus for measurement of transfer functions 

10 in a horizontal direction according to the current invention; 

11 FIG. 3 is a diagram of the coordinate system of used with a 

12 test specimen in the model; 

13 FIG. 4A is a plot of the transfer function magnitude versus 

14 input frequency for the vertical direction test; 

15 FIG. 4B is a plot of the transfer function phase angle 

16 versus input frequency for the vertical direction test; 

17 FIG. 5A is a plot of the transfer function magnitude versus 

18 input frequency for the horizontal direction test; 

19 FIG. 5B is a plot of the transfer function phase angle 

20 versus input frequency for the horizontal direction test; 

21 FIG. 6 is a contour plot of the absolute value of the 

22 dilatational wavenumber on an real -imaginary coordinate system of 

23 the dilatational wavenumbers at 2000 Hz; 
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1 FIG. 7 is a contour plot of the absolute value of the 

2 dilatational wavenumber on an real -imaginary coordinate system of 

3 the dilatational wavenumbers at 5000 Hz; 

4 FIG. 8A is a plot of the real dilatational wavenumber versus 

5 frequency; 

6 FIG. 8B is a plot. of the imaginary dilatational wavenumber 

7 versus frequency; and 

8 FIG. 9A is a plot of the real shear wavenumber versus 

9 frequency; 

10 FIG. 9B is a plot of the imaginary shear wavenumber versus 

11 frequency; 

12 FIG. 10 is a plot of the real and imaginary Young's modulus 

13 versus frequency; 

14 FIG. 11 is a plot of the real and imaginary shear modulus 
15. versus frequency; and 

16 FIG. 12 is a plot of the Poisson ! s ratio versus frequency. 

* 17 

18 DESCRIPTION OF THE PREFERRED EMBODIMENT 

19 The test procedure, consists of vibrating a mass-loaded, 

20 slab-shaped test specimen 10 with a shaker 12 in two different 

, 21 directions, vertical 14A and horizontal 14B, as shown in FIGS. 1 

22 and 2, respectively. It is noted that the load mass 16 attached 

23 to the top of the test specimen 10 must be sufficiently stiffer 

24 than the specimen 10 that it can be modeled as lumped parameter 

25 expression rather than a continuous media system. A typical 
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1 example would be a steel load mass 16 attached above a rubber- 

2 like material test specimen 10. This example results in a ratio 

3 between the two stiffnesses of greater than 100. Lower ratios 

4 result in less accurate estimations. Vibrating the combined 

5 specimen 10 and load mass 16 causes different waveforms to 

6 propagate in the specimen 10. The inverse method developed here 

7 allows for the data from the experiments to be manipulated so 

8 that the complex dilatational and shear wavenumbers. can be 

9 measured for the specimen 10. This test is usually done at 

10 multiple frequencies (swept sine) so any frequency dependencies 

11 can be identified and measured. Input vibration data is 

12 collected from the shaker 12. A sensor 18 is mounted on load 

13 mass 16 and another sensor 20 is mounted on shaker 12 for 

14 collecting transfer function data. In FIG. 1, the test is set up 

15 for monitoring the vertical transfer function. FIG. 2 shows the 

16 test as set up for monitoring the horizontal transfer function. 

17 Sensors 18 and 2 0 should be oriented properly to capture the 

18 motion being measured. Other test configurations using • 

19 directions other than vertical and horizontal are possible; 

20 however, the test setups shown are preferred for ease of set up 

21 and calculation. These sensors 18 can be either accelerometers 

22 that record accelerations, or laser velocimeters that record 

23 velocities. In the swept sine mode, transfer functions of 

24 acceleration divided by acceleration or velocity divided by 

25 velocity are both equal to displacement divided by displacement. 
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1 The time domain data collected from the sensors 18 and 20 are 

2 Fourier transformed into the frequency domain and then recorded 

3 as complex transfer functions, typically using a spectrum 

4 analyzer 22 . 

5 The motion of the test specimen shown in FIGS. 1 and 2 is 

6 governed by the equation: 
7 

8 ^V 2 u + (A + /i)VV«u = p— =- , (1) 

9 

10 where X and ju are the complex Lame constants (N/m 2 ).; 

11 p is the density (kg/m 3 ) ; 

12 t is time (s) ; 

13 • denotes a vector dot product; and 

14 u is the Cartesian coordinate displacement vector of the 

15 material. 

16 The coordinate system of the test configuration is shown in 

17 FIG. 3. Note that using this orientation results in Jb = 0 and a 

18 having a value less than zero. The thickness of the specimen, h, 

19 is. a positive value. Equation (1) is manipulated by writing the 

20 displacement vector u as 
21 



22 



u = 1 



u x (x 9 y 9 z 9 t) 
u y (x 9 y 9 z 9 t) 
[u 2 (x 9 y 9 z 9 t\ 



(2) 
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1 where x is the location along the plate (m) , y is the location 

2 into the plate (m) , and z is the location normal to the plate 

3 (m) , as shown in FIG. 3. The symbol V is the gradient vector 

4 differential operator written in three-dimensional Cartesian 

5 coordinates as 
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^ u . d . d . 

V = — l Y H i v H j_ /q\ 



8 

9 with i x denoting the unit vector in the x-direction, i y denoting 

10 the unit vector in the y-direction, and i z denoting the unit 

11 vector in the z-direction; V 2 is the three-dimensional Laplace 

12 operator operating on vector u as 
13 

14 V 2 u = V 2 u x i x + V 2 u y i y + V 2 u z i z ' (4) 

15 

16 and operating on scalar u as 
17 

■-,2 _ _ 6^ U XVZ <^ U XV7 ^ U Y\)7 

w w 4> 2 ^ 2 ' (5) 

19 

20 and the term V»u is called the divergence and is equal to. 

V7 ^ 2 

dx cy dz ' ■ (6) 
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2 The displacement vector u is written as 
3 



4 
5 



10 

11 
12 
13 
14 
15 

16 

17 
18 
19 

20 



u = + V x , 



(7) 



6 where <j> is a dilatational scalar potential, x denotes a vector 

7 cross product, and y/ is an equivoluminal vector potential 

8 expressed as 
9 



y/ x (x,y,z,t) 
y/ y (x,y,z,i) 



(8) 



The problem is formulated as a two-dimensional system, thus y = 0, 
u y (x,y,z,t) = 0 , and d^)ldy = 0. Expanding equation (7) and breaking 
the displacement vector into its individual nonzero terms yields 



dx. dz 



(9) 



and 



u z (x,z,t)- + 

cz ex 



(10) 
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1 Equations (9) and (10) are next inserted into equation (1) , 

2 which results in 
3 

4 c 2v 2 ^,z,0 = ^^ (ID 

5 

6 and 
7 

c>Vw> = %^ (12) 

9 

10 where equation (11) corresponds to the dilatational component and 

11 equation (12) corresponds to the shear component of the 

12 displacement field. Correspondingly, the constants and c s are 

13 the complex dilatational and. shear wave speeds, respectively, and 

14 are determined by 
15 



i6 c '--f7T (13) 

17 

18 and 
19 

20 c s =..}-■ . (14) 
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1 The relationship of the Lame constants to the Young's and. shear 

2 moduli is shown as 



3 



4 X = ^ (15) 

(l+u)(l-2u) 

5 

6 and 
7 

8 " =G= ^ ' <16) 



10 where E is the complex Young's modulus (N/m ) ,' G is the complex 

11 shear modulus (N/m ), and t> is the Poisson f s ratio of the 

12 material (dimensionless) . 

13 The conditions of infinite length and steady-state response 

14 are now imposed, allowing the scalar and vector potential to be 

15 written as 
16 

1 7 0( x, z, t) = O(z) exp(ita) exp(i cot) (17) 

18 

19 and 
20 

21 y/ y (x,z, t) = ^(z) exp(ifcc) exp(i^) (18) 
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1 where i is the square root of -1, a) is frequency (rad/s) , and k 

2 is wavenumber with respect to the x axis (rad/m) . Inserting 

3 equation (17) into equation (11) yields 
4 



^Q + a 2 ®(z) = 0 , (19) 
dz L 



6 

7 where 
8 



9 <x = ^kl-k 2 , 

10 



(20) 



11 and 
12 

13 *</=—• (21) 

c d . 

14 

15 Inserting equation (18) into equation (12) yields 
16 

17 ^£) +y? 2 V j /(z) = 0 f (22) 

dz 

18 

19 where 
20 

21 J3 = ^kj-k 2 , (23) 
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1 and 
2 

3 k s = — • (24) 
c s 

4 

5 The solution to equation (19) is 
6 

7 0(z) = y4(fc,«)exp(iaz) + exp(-iaz) , (25) 
8 

9 and the solution to equation (22) is 
10 

11 . V(z) = C(k,(o)exp(ij3z) + D(k,co)exp(-i/3z) , (26) 
12 

13 where ,4(fc,a)) , B(k,a>) , C(k,co) , and D(k,co) are wave response 

14 coefficients that are determined below. The displacements can 

15 now be written as functions of the unknown constants using the 

16 expressions in equations (9) and (10) . They are 
17 

u 2 (x,z,t) = U z (k,z,co) exp(ifct) exp(i&rf) 

18 (27) 

= { ia[A(k,eo)exp(iaz)-B(k,(o)exp(-iaz)] + 

ik[C(k,(D)exp(i/3z) + D(k,co)exp(-i/k)\ } exp(ifcc) exp(ia)t) , 

19 

20 and 
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u x (x 9 z,f) = U x (k y z,co) exp(ita) exp(itftf) 

1 (28) 

= { ik[A(k,co)zxp(iaz) + 5(A:,<y)exp(-iaz)]- 

i/?[C(A:,fi))exp(i)fife)-Z)(/r,fl))exp(-iy&)] } exp(i£x) exp(iatf) . 

2 

3 Specific boundary conditions are now needed to individually solve 

4 for the case of vertical and horizontal motion. These are 

5 formulated separately below. 

6 For the case of vertical motion, the base at z = a is 

7 vibrated vertically using a shaker, as shown in FIG. 1. Four 

8 boundary conditions are necessary to formulate this problem. 

9 Because the mass is attached to the material, the tangential 

10 (horizontal) motion at the top of the plate (z = b) is zero and 

11 this equation is written as 
12 

13 u x (x,b,t) = 0 . (29) 

14 

15 The normal stress at the top of the specimen is equal to and 

16 opposite the load created by the mass in the z direction. This 

17 expression is 
18 

19 r Z2 (xM = (A + + x du^M = , M ^u z (xM f (3Q) 

dz dx dt 
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2 

1 where M is mass per unit area (kg/m ) of the attached mass. The 

2 tangential motion at the bottom of the plate (z = a) is zero and 

3 this equation is written as 
4 

5 u x (x,a,t) = Q , . (31) 

6 

7 and the normal motion at the bottom of the plate is prescribed as 

8 a system input. This expression is 
9 

10 u z (x,a,t) = U 0 exp(ia)t) . (32) 

11 

12 Assembling equations (1) - (32) and letting b = 0 yields the 

13 four-by-four system of linear equations that model the system. 

14 They are 
15 

16 Ax = b , (33) 

17 

18 where the entries of equation (33) are 
19 

20 A u =ik , (34) 

21 

22 ^12=^11 > (35) 

23 

24 A l3 =-i/3 , (36) 
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A 14 =~ A \3 • . (37) 



A 2 \ = -a 2 X-2a 2 /u- Ak 2 -\M(o 2 a , (38) 



A 22 = -a 2 l-2a 2 ^-Ak 2 + \Mco 2 a , ■ (39) 



A 23 =-2kfi{t-iMco 2 k , (40) 



A 24 =2kj3^-iMco 2 k , (41) 



A 3l =A n e\p(iaa) , (42) 



^32 = ^1 1 exp(-icoz) , (43) 



,4 3 3 =/t 13 exp(i/k) , (44) 



A ZA =-^13 ex P(-iyfo) , (45) 



A 4l = iaexp(iaa) , (46) 



A 42 =-iaexp(-iaa) , (47) 



A 43 =ikexp(i/3a) , (48) 
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1 /I44 =i£exp(-i/fc) , (4 9) 

2 

3 x n =A(k,o)) , (50) 

4 

5 x 2 i=5(£,<y) » (51) 

6 

7 * 31 =C(*,6>) , (52) 

8 

9 x 41 , .-.'.* (53) 

10 

11 ^11=0 , (54) 

12 

13 b 2l =0 , (55) 

14 

15 * 31 =0 , (56) 

16 and 
17 

18 &41=^0 • (57) 

19 

20 Using equations (34) - (57) the solution to the constants A(Jc,co) , 

21 B{k,a>) , C{k,ai) , and D(k,co) can be calculated at each specific 

22 wavenumber and frequency using 
23 

24 x = A _1 b . (58) 



17 



1 Noting that for vertical motion, k = 0, and using the 

2 coefficients from equation (58) , the transfer function between 

3 the vertical base displacement and the vertical mass displacement 

4 can be written as 



5 



Ji(a>) = 



l u 2 (0M l 



v cos(k d h)- 



— \k d sm(k d h) 



(59) 



7 

8 where Ti(co) or R\(o)) correspond to the data from the vertical 

9 motion experiment. 

10 The next step is to solve the inverse problem for vertical 

11 motion. This involves using the experimental data and equation 

12 (59) to estimate the dilatational wavenumber. Equation (59) can 

13 be rewritten as 
14 

15 f(k d ) = 0 = cos(k d h) - \^k d sm{k d h) - R x (60) 

16 

17 where the problem now becomes finding the zeros of the right-hand 

18 side of equation (60), or, in the presence of actual- data that 

19 contains noise, finding the relative minima of the right-hand 

20 side of equation (60) and determining which relative minimum 

21 corresponds to dilatational wave propagation and which relative 

22 minima are extraneous. Because equation (60) has a number of 

23 relative minima, zero. finding algorithms are not applied to this 
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1 function, as they typically do not find all of the minima 

2 locations and are highly dependent on initial starting locations. 

3 The best method to find all of the minima locations is by 

4 plotting the absolute value of the right-hand side of equation 

5 (60) as a surface with the real part of dilatational wavenumber 

6 k d on one axis and the imaginary part of k d on the other axis. 

7 In order to do this, the user should start at a low frequency 

8 where the aliasing minimum has not yet appeared. In the specific 

9 example shown herein, this is below 3850 Hz for the dilatational 
. 10 wave and below 1550 Hz for the shear wave. At these lower 

11 frequencies, the minimum furthest to the left will correspond to 

12 dilatational wave propagation. As the frequency increases, 

13 extraneous minima will appear to the left of the minimum that 

14 corresponds to dilatational wave propagation, however, the wave 

15 propagation minimum will always be close to the previous test 

16 frequency wave propagation minimum provided that the frequency 

17 increments are relatively small. At a resolution of 0.5 rad/m 

18 for the materials in the example herein, this requires a 

19 frequency increment of 37.3 Hz for the dilatational measurement 

. 20 and 14.4 Hz for the shear measurement. Different test specimens 

21 and top masses require different increments. Additionally, the 

22 real part of the wavenumber is a monotonically increasing 

23 function with respect to frequency, so at each increase in 

24 frequency, the new wavenumber to be estimated has to be greater 

25 than the old wavenumber that was previously estimated. This 
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1 process is further illustrated as related to the discussion 

2 concerning FIG. 6 and FIG. 7 below. 

3 For the case of horizontal motion, the base at z = a is 

4 vibrated horizontally using a shaker, as shown in FIG. 2. Four 

5 boundary conditions are necessary to formulate this problem. 

6 Because. the mass is attached to the material, the shear 

7 (tangential) stress at the top of the plate is equal to .opposite 

8 the load created by the mass in the x direction. This expression 

9 is 



10 

11 Tzx( x M = P 

12 



du x {xfrt) t du 2 {xfrt) 
dz 3c 



= _ M ctu x (x M t (61) 



a 1 



13 where M is mass per unit area (kg/m 2 ) of the attached mass. The 

14 normal motion at the top of the plate (z = b) is zero and this 

15 equation is written as 
16 

17 u z (x,b,t) = 0 . . (62) 

18 

19 The tangential motion at the bottom of the plate (z = a) is 

20 prescribed as a system input and this equation is written as 
21 

22 u x (x,a,t) = V 0 exp(ici)t) , (63) 
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1 and the normal motion at the bottom of the plate is zero. . This 

2 expression is 
3 

4 u z (x,a,t) = 0 . (64) 

5 

6 Assembling equations (1) - (28) and (62) - (64) and letting 

7 b = 0 yields the four-by- four system of linear equations that 

8 model the system. They are 
9 

10 Ax = b , (65) 

11 

12 where the entries of equation (61) are 
13 

14 A n =-2fjka-\co 2 Mk , (66) 

15 

16 A l2 =2iAa-'\a> 2 Mk , (67). 

17 

18 A u =fiip 2 -vk 2 + i6) 2 Mp , (68) 

19 

20 A l4 =^ 2 -^k 2 -id) 2 Mfi , ■ (69) 

21 

22 A 2i =ia , (70) 

23 

24 A 22 =-A 2[ , (71). 
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1 

2 A 23 =ik , (72) 
3 

4 ^24=^23 . (73) 
5 

6 A 3l = A 23 exp(iaa) , (74) 
7 

8 A 32 =A 23 exp(-iaa) , (75) 
9 

10 ^33=-i>3exp(i/&) , (76) 

11 

12 A 34 =i/3exp(-ij3a) ,' (77) 

13 

' 14 ^41 = -^21 exp(iaa) , (78) 

15 

16 A 42 =-A 2l exp(-iaa) , (79) 

17 

• 18 A 43 = A 23 exp(ij3a) , (80) 

19 

20 A 44 =A 23 Qxp(-i/h) , (81) 

21 

22 

23 x u =A(k,a>) , (82) 

24 

25 x 2 l , (83) 
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1 

2 x 3l =C(k,G)) , (84) 
3 

4 x 4l = D(k,a)) , (85) 

5 
6 

7 fcn =0 , (86) 
8 

9 &2i=0 , (87) 
10 

11 & 31 = F 0 ' < 88 > 

12 and 

13 

14 b 4l =0 . (89) 
15 

16 Using equations (67) - (89) the solution to the constants A(k,co) , 

17 B(k,a>), C(k,a>) , and LHJk,a>) can be calculated at each specific 

18 wavenumber and frequency using 
19 

20 x = A _1 b . (90) 
21 

22 Noting that for horizontal motion, k = 0, and using the 

23 coefficients from equation (90) , the transfer function between 

24 the horizontal base displacement and the horizontal mass 

25 displacement can be written as 

23 



8 



r2(g) ,^-,"*<ftM).. 



R ^ V ° cos(M)- 



(91) 



^sin^A) 



2 

3 where T 2 ((o) or correspond to the data from the horizontal 

4 motion experiment. 

5 The next step is to solve the inverse problem for horizontal 

6 motion. This involves using the data and equation (91) to 

7 estimate the shear wavenumber. Equation (91) can be rewritten as 



f(k s ) = Q = cos(k s h)- 



— \k s sm(k s h)-R 2 . (92) 



10 

11 It is noted that this equation is identical, except for the 

12 subscripts, to equation (60) . The shear wavenumber is estimated 

13 using the same procedure that was used to estimate the 

14 dilatational wavenumber above. 

15 The material properties can be determined from the 

16 wavenumbers. First, the dilatational and shear wavespeeds are 

17 determined using 
18 

!9 c d=T~ 03) 

k d 

20 

21 and 
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1 C s =^- , (94) 

2 

3 respectively. The Lame constants are calculated using equations 

4 (13) and (14) written as 
5 

6 M = pCs (95) 

7 

8 and 
9 

10 X = pc 2 d -2pc 2 . (96) 

11 Poisson's ratio is then calculated using 
12 

13 v= . (97) 

2Cu+A) 

H 

15 Young's modulus can be calculated with 
16 



17 £ _ 2//(2//+3A) 
' . " 2(//+A) 



(98) 



18 

19 and the shear modulus can be determined using 
20 

21 G = m • (99) 

22 The above measurement method can be simulated by means of a 

23 numerical example. Soft rubber-like material properties of the 

25 



1 test specimen are used in this simulation. The material has a 

2 Young's modulus E of [(le8-i2e7) + (5e3/-i3e2/)].N/m 2 where f is 

3 frequency in Hz, Poisson's ratio u is equal to 0.40 

4 (dimensionless) , density p is equal to 1200 kg/m 3 , and a 

5 thickness h of 0.1 m. The top mass is a 0.0254 m (1 inch) steel 

6 plate that has a mass per unit area value M of 199 kg/m 2 . FIG. 4 

7 is a plot of the transfer function of the system for vertical 

8 motion and corresponds to equation (59) . FIG. 5 is a plot of the 

9 transfer function of the system for horizontal motion and 

10 corresponds to equation (91). In FIGS. 4 and 5, the top plot is 

11 the magnitude in decibels and the bottom plot is the phase angle 

12 in degrees. 

13 FIG. 6 is a contour plot of the absolute value of equation 

14 (60). expressed in decibels versus real dilatational wavenumber on 

15 the x axis and imaginary dilatational wavenumber on the y axis at 

16 2000 Hz.. The estimated dilatational wavenumber, read directly 

17 from the plot at the location the minimum value appears and 

18 marked with an arrow, is 27.89 + 2.61i rad/m. The actual value 

19 of the dilatational wavenumber is 27.99 + 2.60i rad/m, which is 

20 slightly different from the estimated value due to the surface 

21 discretization of equation (60) . FIG. 7 is a contour plot of 

22 equation (60) at 5000 Hz . At this frequency, an extraneous 

23 minimum appears on the left-hand side of the plot. However, 

24 because the real part of the wavenumber must be increasing with 
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1 increasing frequency, the minimum corresponding to dilatational 

2 wave propagation is located at the arrow marked spot and is equal 

3 to 65.86 + 5.60i rad/m, as compared to an actual value of 65.77 + 

4 5.62i rad/m. Again, the difference between the two values can be 

5 attributed to the discretization of the surface. 

6 FIG. 8 is plot of actual (solid line) and estimated (x 

7 symbols) dilatational wavenumber versus frequency. FIG.. 9 is 

8 plot of actual (solid line) and estimated (+ symbols) shear 

9 wavenumber versus frequency. In FIGS. 8 and 9, the top plot is 

10 the real part of the wavenumber and the bottom part is the 

11 imaginary part of the wavenumber. FIG. 10 is a plot of actual 

12 (solid line) and estimated (real part - x symbols, imaginary part 

13 - o symbols) Young's modulus versus frequency. FIG.. 11 is a plot 

14 of actual (solid line) and estimated (real part - x symbols, 

15 imaginary part - o symbols) shear modulus versus frequency. In 

16 FIGS. 10 and 11, the imaginary part of the modulus all have a 

17 negative sign but are depicted with positive signs for plotting 

18 purposes. FIG. 12 is a plot of actual (solid line) and estimated 

19 (square symbols) of the real part of Poisson's ratio versus 

20 frequency. Because the numerical example is formulated using a 

21 Poisson's ratio that is strictly real, no imaginary component is 

22 shown in this plot. Imaginary values of Poisson's ratio are 

23 possible and have been shown to theoretically exist. 

24 This invention gives the ability to estimate complex 

25 dilatational and shear wavespeeds of a material that is slab- 
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shaped and subjected to compressive forces. It also allows 
estimation of complex Lame constants of a material that is slab- 
shaped and subjected to compressive forces: Complex Young's and 
shear moduli of a material that is slab- shaped and subjected to 
compressive forces can be estimated using this invention. The 
invention also allows estimation of the complex Poisson' s ratio 
of a material that is slab- shaped and subjected to compressive - 
forces.. The advantage of this patent is that it does not require 
a testing configuration that has to be pressurized. 

Obviously many modifications and variations of the present 
invention may become apparent in light of the above teachings. 
In light of the above, it is therefore understood that within the 
scope of the. appended claims, the invention may be practiced 
otherwise than as specifically described. 
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